% 创建时间：2025-06-21 14:02:32
% 作者：兪馳
clc;
clear;
A = [-2, 1, -3; -1, 0, -2; 0, -1, -2];
B = [2, 2, -2; -1, 5, -2; -1, 1, 2];
C = [5, -4, 4; 1, 0, 4; 1, -1, 5];
X = are(A, B, C)
norm(A' * X + X * A - X * B * X + C)